Pyramidal Atoms: Berylliumlike Hollow States. 
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Based on the idea that four excited electrons arrange themselves around the nucleus in the cor- 
ners of a pyramid in order to minimize their mutual repulsion, we present an analytical model of 
quadruply excited states. The model shows excellent comparison with ab initio results and provides 
a clear physical picture of the intrinsic motion of the four electrons. The model is used to predict 
configuration-mixing fractions and spectra of these highly correlated states. 

PACS numbers: 31.10.+z,31.25.-v 



Describing electron correlations in terms of simple 
physical pictures requires an understanding of the inter- 
electronic interaction which goes beyond numerical di- 
agonalization of the Hamiltonian of the system. The 
study of electron-electron correlation in multiply excited 
states has contributed a lot to this understanding and 
in the eighties group theoretical P, molecular adia- 
batic 01 , and hyperspherical Q approaches led to an un- 
derstanding of doubly excited states in terms of approxi- 
mate quantum numbers associated with angular and ra- 
dial correlations. In the nineties advances in light sources 
stimulated research on triply excited hollow states in Li 
(see Ref. 01 for a recent review). In particular, theoret- 
ical models referred to as the three-electron ionization 
ladder 0, Q, the normal mode analysis 0, 0, and 
the symmetric rotor (ill fl^ IT^ confirmed the appealing 
physical picture for the lowest-lying triply excited states 
that the three electrons be arranged in an equilateral 
triangle coplanar with the nucleus situated in the cen- 
ter, and identified the magnetic quantum number with 
respect to the body-fixed z axis as approximately con- 
served. 

Current interest in multiply excited states is strongly 
stimulated by further development of synchrotron radi- 
ation sources and free-electron lasers providing unprece- 
dented brightness in the uv and xuv and therefore read- 
ily access to the spectral region of interest. Success- 
ful access to hollow lithium states j3| makes studies of 
quadruply excited states a natural next step in this de- 
velopment. Sofar there are no experiments and only very 
few theoretical studies on these systems. One work [lj| 
considered the lowest-lying quadruply excited state 
in Be by multiconfigurational Hartree-Fock calculations 
while two other works (iH presented a classification 
scheme based on symmetry considerations. Hence, in 
view of experimental advances, it seems timely to aim 
at a better understanding in the form of physical im- 
ages, including geometry, symmetry considerations and 
electronic motion for highly correlated four-electron sys- 
tems. Here we shall concentrate on intra-shell states with 
the electrons confined to the same shell. These states 
are interesting since the dominant configurations of the 
lowest-lying quadruply excited beryUiumlike states are 



intra-shell states. 
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Figure 1: Thresholds and the lowest quadruply excited P'' 
state in Be. The insert to the right shows, in the body-fixed 
frame, the physical picture on which our model is based. The 
nucleus is located in the origin and the circles denote the 
electrons. 



Figure shows the thresholds in Be, the lowest '^P'^ 
quadruply excited resonance state and an insert of the 
physical picture on which our model is based. We note 
that while two and three electrons are confined to a line 
and a plane, respectively, the extent of four randomly 
placed electrons is three dimensional. Furthermore, the 
electron-electron correlations are expected to be stronger 
when more electrons interact. Therefore, research on 
these states opens new possibilities to test theoretical ap- 
proximations for correlated electron motion. 

We consider the Coulombic system of four electrons 
moving in the field of a point-like nucleus of charge Z 
corresponding to the nonrelativistic Hamiltonian 
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where fi is the reduced electron mass and e is the elec- 
tron charge. The electron-electron interaction compli- 
cates the problem since the potential is not diagonal in a 
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basis of antisymmetrized products of single-electron con- 
figurations. Consequently, single-electron quantum num- 
bers are not conserved and the eigenfunctions are rep- 
resented as superpositions of configurations. Our model 
accounts for the electron-electron interaction by assum- 
ing the wave function be a product of carefully chosen 
one-electron states. Furthermore, the ansatz for the wave 
function must be supplemented to take into account basic 
quantum symmetries (parity, spin, etc.) 

The Coulomb repulsion between the four electrons is 
minimized when they are placed in the vertices of a reg- 
ular tetrahedron (see Fig. In this configuration the 
interaction between any electron and the remaining three 
is represented by a repulsive electric field along the direc- 
tion, ^, of the electron considered with respect to the nu- 
cleus. Hence, the Hamiltonian in Eq. 1^ can be approx- 
imated by a sum of hydrogenic Hamiltonians perturbed 
by weak external electric fields causing only intra-shell 
transitions. The eigenstates of the hydrogenic Hamilto- 
nian perturbed by a weak external electric field are the 
polarized Stark states, \nkm), described by the princi- 
pal, parabolic and azimuthal quantum numbers denoted 
n,k and m, respectively The wave function there- 

fore splits up into a product of one-electron states which 
is the simplest possible way to represent a many-particle 
wave function. Thus, the basic idea in the model is that 
the inter-electronic interaction tends to stabilize the elec- 
trons in individual Stark states (see Refs. for two 
and for three electron cases) . The parabolic quantum 
number, k, may be associated with the electric dipole mo- 
ment of the atom through the relation pz = ^enk [l^ . 
Hence, the largest energy shift and thus the lowest-lying 
state for each n is obtained when the Stark states are 
maximally polarized along \n,k = n — l,m = 0)^^. 
The ansatz for the wave function is therefore 

\Tr) = \nn ~ 10) ^,\nn - 10) ^,\nn- 10) ^,\nn- 10)^, , (2) 

with the direction of the jth electron with respect to 
the nucleus. Quantum numbers in kets without a sub- 
script refer to quantization along the z axis so the elec- 
tron on this axis is represented by the state \nn — 10). 
The other three states, \nn — 10)^^, q = 1,2, 3, are related 
to it via rotations \nn — 10)^^ — Il[q^, $, O) \nn — 10) = 

e-"}^^'e-"^^y\nn - 10) , with $ = 109.5° the inter- 
electronic angle as seen from the nucleus (see Fig.Q. 

The wave function \Tfi) is anisotropic since the elec- 
trons are oriented. As a consequence of this break of 
rotational symmetry, |Ti^) is not an eigenstate of and 
1/2 . Contrary, the rotational average over Euler angles 
uj = (a,/3,7) 



I^LMAf,.^^ / da d7 / dp sin (3Dijj^^{u;r\T n)^ , 
Jo Jo Jo 

(3) 

is so [121 . Here n denotes the parity and Mj, is the 
magnetic quantum number with respect to the body- 



fixed z axis. In molecular theory this is an important 
quantum number specifying the rotational state of the 
molecule. Since the Hamiltonian is spin-independent the 
total wave function is simply a product of the spatial 



function, |T^*'*'^''), and a spin function, Ixm! )• ^he 
spin function specifies the total spin, S, of the four elec- 
trons and its projection, Ms, which are both conserved 
quantities. The indices, 5*12 and S'34, refer to the interme- 
diate spins arising from coupHng the spin of the first and 
second and the third and fourth electron, respectively. 
The total wave function is obtained by antisymmetrizing 
the product state \T'^'^^"')\xMf "^'^) , 
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where A is the four-particle antisymmetrization opera- 
tor. The final expression for the four-electron state is 



LMMiTT.SMs 



N being a normalization constant. The possible values 
of intermediate spins, iS'12 and iS'34, are determined by 
S — S12 + S'34 added as angular momentum operators. 

Our model provides an analytical expression for the 
wave function describing quadruply excited intra-shell 
states of the four-electron atom and gives a simple 
physical image of the atomic states. To test the 
quality of the model, configuration-mixing fractions 
are compared with ah initio results. We expand 
the ansatz in the spherical basis ^ii^^^'^i'^^^^'^') = 
Szpm, c{lfiThs, LMMjTT, SMs) IkfJ'imsi), {p = 

(Pl,P2,P3,P4),P = {l,fi,ms)) with \-q,LMMi^^SAh^ ^^^^^ 

normalized such that X^F/ime \'^i^P'''^s, LAIMjtt, SMs)\'^ — 
1. The configuration-mixing fractions specifying the I 
mixing between several multiplet states are defined as 



with 



pLMnrS 
la 



J2 J2 \cilfifhs,LMMnr, SMs)\^ , (5) 



where Ic = {hhhlijc denotes the /-configuration 
regardless of permutation, e.g., {sppp)c — 
{sppp,pspp,ppsp,ppps}. For n = 2 the pyramidal 
states are shown in Figl^b). The mixing fractions 
for and ^D*' are 0.9 s^p^ and 0.1 p^. The other 
terms have only a single configuration (sp'^). Table 
ID shows a comparison between the mixing fractions 
predicted within the model and from multiconfigura- 
tional Hartree-Fock calculations The model not 

only predicts the trend in the mixing fractions, but 
actually compares very well with the configurational 
calculations. The convincing agreement in the pre- 
dictions shows that our four-electron model accounts 
for the strong electron-electron correlations within the 
quadruply exited intra-shell states. Table shows 
configuration-mixing fractions for all terms with L < 4 
within the n = 3 shell for which the pyramidal shape 
is accessible [iH ll^ . The mixing fractions are seen 
to be strongly spin dependent. As a consequence of 
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Table I: Configuration-mixing fractions in Be for the most important intra-shell configuration states of symmetry. SR denotes 
predictions within the present symmetric rotor model, and MCHF denotes multiconfigurational Hartree-Fock calculations 0|. 



n=3 n=4 n=5 n=6 



Configuration 


SR 


MCHF 


SR 


MCHF 


SR 


MCHF 


SR 


MCHF 


sp^ 


0.94 


0.90 


0.81 


0.76 


0.67 


0.66 


0.55 


0.59 


spd^ 


0.063 


0.084 


0.15 


0.18 


0.20 


0.21 


0.23 


0.23 


sd^f 






0.013 


0.0081 


0.039 


0.040 


0.068 


0.053 


p^df 






0.022 


0.032 


0.057 


0.058 


0.091 


0.078 


pdf^ 






0.0014 


0.0025 


0.0085 


0.0081 


0.021 


0.017 


d^f 






0.0022 


0.0025 


0.0097 


0.0064 


0.021 


0.012 



Pauli statistics some configurations are even forbidden 
for accessible terms. Furthermore, there is a notable 
difference between the mixing fractions for quadr uply 
excited states and those for triply excited states |l2l |: 
Contrary to triply excited states, the mixing fractions 
for quadruply excited intra-shell states do not depend 
on Ml. This refiects the three dimensional symmetry 
of the regular tetrahedron (spherical top) compared to 
two dimensional symmetry of an equilateral triangle 
(symmetric top). 

The energy of a state within the symmetric rotor model 
is given by (vp^*^^'^^ ItjI vj/^*^^-^), i.e.. 
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with rij = |7\ — rj|. Here En is the Bohr energy and the 
last term is the electron-electron interaction energy. The 
equality follows from the fact that the Stark states are 
eigenstates of the single-particle part of the Hamiltonian 
of Eq. (0. 

Figure|2shows the energies predicted within our model 
for the pyramidal states within n — 2 and for L < 4 
within the n = 3 shell. The energy increases with in- 
creasing L. Furthermore, we see the tendency that the 
energy increases with increasing spin, S, contrary to the 
behavior in ground states configurations where the term 
with the largest possible value of S for a given L has low- 
est energy (Hund's rules). Finally, the energy levels split 
up for different parity since different parity requires dif- 
ferent configurations. We note that if the system is con- 
sidered from a molecular viewpoint the M/-independence 
of the configuration-mixing fractions is expected since 
the rotational energy of the spherical top molecule reads 
E — jjL{L + 1), with / — X^i^^^f being the moment 
of inertia of the molecule. However, the L scaHng does 
not follow this relation and it is clear from Eq. © that 
rotationallike structure in excited atomic states is not as- 
sociated with the kinetic energy term as in a molecule, 
but is exclusively due to the electron-electron interaction. 



The term, ( * 
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R 



the nuclear charge, Z, due to the scaling of the matrix 
element of — . Contrary, the spherical top model pre- 

diets a Z"^ scaling of the energies due to the scaHng of 
the moment of inertia. 

In summary, we constructed an analytical wave func- 
tion for quadruply excited states corresponding to the 
case where the four excited electrons are in the corners 
of a pyramid. By comparison with ah initio calculations, 
we confirmed the accuracy of the model and, hence, the 
appeahng physical image of the correlated motion of the 
four electrons. We predicted configuration mixing frac- 
tions for a number of states and generally found that a 
labeling of states in terms of single configurations is insuf- 
ficient. Finally, we calculated the spectrum of pyramidal 
states within the second and third principal shell and 
compared with the expectations obtained when viewing 
the system as a molecularlike spherical rotor. 

We propose to address pyramidal states by dipole 
allowed light absorption transferring a fraction of the 
Be ground state population to the singly excited 
Be(lsls2s2p^P°) state by resonant laser excitation (5.277 
eV). A synchrotron or a free electron laser source could 
then subsequently scan through the region of quadruply 
excited pyramidal states of ^D*' symmetry. A similar two- 
color technique was successfully applied on triply excited 
states in Li [201 . 
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Table II: Configuration-mixing fractions for pyramidal states 
in Be for the n = 3 shell with L < 4. 
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Figure 2: Energy levels for pyramidal intra-shell states in Be. 
(a) n = 3 states with L < 4. (b) n = 2 states. 
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